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Abstract 

We compute the parametric correlation function of the conductance peaks in 
chaotic and weakly disordered quantum dots in the Coulomb blockade regime 
and demonstrate its universality upon an appropriate scaling of the param- 
eter. For a symmetric dot we show that this correlation function is affected 
by breaking time-reversal symmetry but is independent of the details of the 
channels in the external leads. We derive a new scaling which depends on 
the eigenfunctions alone and can be extracted directly from the conductance 
peak heights. Our results are in excellent agreement with model simulations 
of a disordered quantum dot. 
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Quantum dots |T],|| are semiconductor heterostructures that confine about 100 electrons 
by an electrostatic potential to isolated islands with a typical size of less than a micron. 
The dots can be weakly coupled via tunnel barriers to external leads in order to study their 
transport properties. For sufficiently low temperatures the conductance of the dot exhibits 
equally spaced peaks with increasing gate voltage, where each successive peak corresponds 
to a tunneling of a single electron into the dot. This occurs when the increase in the 
Fermi energy in the leads matches the energy required to charge the dot by one additional 
electron. The suppression of tunneling between the peaks by Coulomb repulsion is known as 
Coulomb blockade. A striking feature of these resonances is the irregular dependence of their 
amplitudes on controllable parameters, such as the shape of the dot or an external magnetic 
field. These fluctuations have recently been accounted for by a statistical theory based on 
the assumption that the Hamiltonian of the dot can be described by random-matrix theory 
(RMT) due to irregularities in the confining potential which give rise to a chaotic classical 
dynamics. RMT description is suitable also for dots with weak impurity disorder since the 
same conductance distribution is obtained when the dot is modelled by a random potential 

i- 

Recent experiments have been probing the conductance of quantum dots as an external 
parameter is varied. In this paper we study the correlation between conductance peak am- 
plitudes which belong to different values of the parameter. Very little is known about this 
parametric correlator in the regime of isolated resonances, in contrast with the overlapping 
resonance regime characteristic of open dots where correlations versus energy and magnetic 
field were investigated both theoretically (for a large number of channels) and experi- 
mentally . By casting the peak correlator in the framework of the Gaussian random-matrix 
process (CP) ||^, we demonstrate its universality upon an appropriate scaling of the param- 
eter and compute its universal form. Our main results for the conductance correlator are 
the approximate expression (|^) and its exact short- distance behavior ([1^. We also derive 



an alternative parameter scaling which is extracted directly from the measured conductance 
peaks according to Eq. (|T5|). Its universal ratio to the usual level velocity scaling (||) is 



given by Eqs. (0) and ([T2D for the orthogonal and unitary cases, respectively. 

We focus on the temperature regime F <C kT < A, typical of experiments where 
A is the mean resonance spacing. In this regime only one quasi-bound state (usually the 
ground state) of the electron gas in the dot contributes to each conductance peak |10|,|TT|, 
whose width is ~ kT. The peak amplitude that corresponds to a level A is given by [|r^ 

r ^ ^aF-a ... 
^ h 2kT r';, + ' ^ ^ 

where r''jf^ is the total decay width from level A into the left (right) lead. Each lead can 
support several open channels so that F^^*^^ = J2c'^'^}x\ where F^.^^^ is the partial decay width 
into channel c. A suitable formalism that relates the level width to the dot's eigenfunction 
is the i?-matrix theory [|12[, according to which F^a =| 7cA P where 



7cA 



{h\,Pjmy^^ J dS<5l^x (2) 



is the partial width amplitude, kc is the channel momentum and Pc is a penetration factor 
to tunnel through the barrier. Here \E'a is the dot eigenfunction, $c is the transverse wave- 
function in the lead corresponding to an open channel c, and the integral extends over the 
cross-section of the lead at its end attached to the dot. 

The connection with RMT is made by assuming that the statistical properties of 
the dot's Hamiltonian are well described by a x random matrix H taken from the 
appropriate Gaussian ensemble (GE). The dot's eigenfunction \1/a is related to the A-th 
eigenvector of H, -ipx^, via the expansion '^x{r) = Y^^=iip\^p^{r), where p^{r) are a set of 
solutions of Schrodinger's equation at energy Ex inside the dot. Omitting the subscript A, 
the partial width amplitude (|]) can be written as a scalar product 7c = I]/^ = (0c | V') , 

where (p^^ = {h'^kcPc/'mY^'^ J dS^cP*^- The correlations between the 7c define the channel 
correlation matrix M, and are given by the scalar products of the channel vectors M^c' = 
7*7c' = jf{4>c I (pc')- The channels can have any degree of correlation among them and 
different decay widths into them, hence 0c are in general not orthogonal and have different 
norms. 



In order to discuss the correlation between conductances at different shapes or external 
fields we use the framework of the GP, which generalizes Dyson's GE to describe statistical 
properties of systems that depend on a parameter P,p^. A GP is a set of random N x N 
matrices H{x) whose elements are distributed at each x according to the appropriate GE 
with a prescribed correlation among elements at different values of x: 



HxJx) = 



HUx)H,^{x') = ^f{x - x')gi'l^ , (3) 

where g'xa~i2 = ^X/i^au + ^Xv^ufi for the Gaussian orthogonal process (GOP) and g^\~^} = 
'^^Xu^cTfi for the Gaussian unitary process (GUP). Assuming that the leading order of / is 



/ ^ 1 — k(x — x'Y, we have shown p that k = P^jf{dEx/dxy/A'^ which suggests scaling 
X by the RMS of the level velocity as originally discussed in 0: 

'{dEx/dxY/A'^] ^'^ X = VDx . (4) 



X X 



After the scaling (^) / ^ 1 — (3^{x — x'Y/N becomes independent of the non-universal 
quantity k. The parametric correlator for any observable, being determined by A^(l — /), is 
then a universal function of x — x' . This applies not only to spectral correlators but 
also to those involving the eigenfunctions. In particular, the conductance peak correlator 

G{x)~G 



_ ^ G{x)G{x') , G{x) = (5) 

becomes universal upon the scaling (^) for all dots characterized by the same channel cor- 
relation matrix M. 

The case of a left-right symmetric dot is particularly simple since = F*" = F and the 



conductance peak correlator reduces to the width correlator F(x)F(x'). We note that the 
correlation matrix M^c' is Hermitean and positive definite, and hence can be transformed into 
a diagonal form which defines a set of orthonormal eigenchannel vectors (pc- Expressing the 
width in terms of the partial widths of the normalized eigenchannels fc =| (0c | i^) and 
using the property that the cross-channel correlator is smaller by than the autochannel 



correlator, we obtain in the limit of large N [14 



CG{x-x') = t,{x)t^{x') , (6) 

where Ti is the partial width of an arbitrary normahzed channeL Remarkably, the peak 
correlator for a symmetric dot is not only universal but is also independent of the details of 
the channels in the external leads, including their number, the rate of tunneling into each and 
the correlations among them. We computed this function for both the orthogonal and the 
unitary cases from simulations of the simple GP defined by [|T^ H{x) = Hi cosx + H2 sinx 
where H12 are independent GOE (GUE) matrices. The results are presented in Fig. ^ We 
found that cq is fitted very well by a Lorentzian in the orthogonal case (/3 = 1) and by a 
squared Lorentzian in the unitary case (/5 = 2): 



cg{x - x') 



(7) 



1 + (x — x'Y jo^p 

where a\ = 0.48 ± 0.04 and 0:2 = 0.64 ± 0.04 with a P^r degree of freedom of ~ 10~^. 

We tested the GP prediction by studying the peak correlator in a disordered dot modelled 
by a two-dimensional Anderson Hamiltonian. We used a 27 x 27 lattice with a cylindrical 
geometry and introduced a parametric dependence by adding a step potential of strength 
X along the symmetry axis. The left and right leads are represented by arrays A^^"^^ of 
mi X 1712 lattice points with total widths given by F''^''^ = Y,r,€A^<~^) "^I I ^(^«) Fig. |I| 

we show the peak correlator for 1x1 and 4x4 leads, corresponding to leads with a single 
channel and with many correlated channels, respectively. We took Vi = 1 but verified that 
different choices do not change the results, which are in excellent agreement with the GOP 
prediction. In the inset we show F(a;)/r for a typical member of the GOP and for one 
realization of the site energies in the Anderson model, using a single eigenfunction without 
any statistical averaging. As expected, this quantity exhibits irregular oscillations with a 
period comparable to the decorrelation distance along x. 

In order to break time-reversal symmetry, we applied a constant magnetic field along the 
symmetry axis, tuned such that the cylinder encloses a fiux of 1/4 fiux unit. The parametric 
dependence was introduced either by a step potential as before or by closing the cylinder 
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into a torus and applying a varying magnetic field perpendicular to the constant one. The 
peak correlators in both cases are also plotted in Fig. |1|, and agree very well with the GUP 
prediction. 

Although in principle all parametric correlators become universal upon the scaling 
in practice this is not always experimentally feasible. The energy levels probed in a quantum 
dot are usually not the excited states for a fixed number of electrons, but the ground states 
for different numbers of electrons. Since the spacing between the conductance peaks is 
dominated by the charging energy e^/C which is much larger than the mean spacing A, the 
scaling factor in (^) is difficult to measure. It is therefore important to derive a scaling that 
can be extracted directly from the conductance peak heights. For that purpose we define 
■j\{x) = (0 I ip\{x)) and consider its derivative 

dx I'-x x'- x^^ Ex-E^ "" ^ ' 

given by first-order perturbation theory and where unprimed quantities are evaluated at 
X. The RMS of d'jx/dx is then calculated in two steps. First, we perform the average 
over H{x') at fixed H{x), employing the conditional probability distribution |^ for which 



Hf^x{x')H*^{x') = 6f^,^a'^{l — P)/P- Next we average over H{x), taking advantage of the 
factorization of P [H{x)] into a product of the eigenvalue and eigenfunction probability 
distributions. We obtain the following expression for the non-universal quantity k in the 
expansion of / preceding (Hf): 



1 , d'jx 



where the constant Z = (aV/5A^)E^^A(^A - E^)-^ = (4/^7^^) /o°° de[l - l2(e)]/e^ and F2(e) 
is the two-point cluster function. Eq. (Pf) suggests a new scaling: 



1/2 



^-'^r=\^^\^\y\lX?j X^VRx. (10) 

With this scaling, / 1 — (xj. — x'^)^/ZN, and all correlators become universal functions 
of Xr — x'^- The scaling factor in ( p!0D can be interpreted as the RMS of the eigenfunction 
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rotation rate in analogy with the energy level velocity in the scaling (^). The quantity Z 
above diverges for the GOP but is finite for the GUP Z = 2/3, hence the scaling (p!0| ) is 
well-defined only in the unitary case. In this case (GUP) the ratio between the rotation rate 
and level velocity scaling factors is a universal constant, obtained from comparing with 
the expression for k preceding (^) 

R/D = ttVS (GUP) . (11) 

The divergence of the scaling factor in (p]0| ) for the GOP notwithstanding, it is still pos- 
sible to use this scaling if the derivative d'jx/dx is replaced by A^x/^x. This regularization 
translates into a small-spacing cutoff 6 in Z such that \ Ex — \ /A < 6 are excluded from 
the sum, resulting in a logarithmic divergence Z oc log 5. For a small but finite Ax we can 



then deduce the ratio between the two scaling factors for the GOP fT^ : 

R/D = log(Ax) + Const. (GOP) . (12) 
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Relation ( |T^ ) is a universal function of Ax, as we demonstrate in the inset of Fig. and as 
is also confirmed in simulations of the Anderson model. 

We now come back to the conductance peak correlator (^ and extract its short-distance 
behavior. We consider the quantity (AF^)^ = [FA(a;') — Fa(x)]^ and calculate its average in 
first-order perturbation theory by expanding ipx{x') in iIj^{x). Averaging first over H{x') we 
obtain 

(ATIF ^ ^(1 - n E [^A - E.r'TxT, , (13) 



where the remaining average is over H{x). Using the GE relation F^F^ = 
[(3N/2{N - 1)] (Ff - Tx^) and the scaling (|10D we find 

R 

cg{x - x') ^ 1 - hp— I X - x' p , (14) 

where hp = {3 and R/D is given by (0) for the GUP and by (|g) for the GOP. When 
compared with the leading-order behavior of the squared Lorentzian (^ for the GUP case 



we find a2 = VS/n ~ 0.55. The discrepancy with the value quoted below (|^ indicates 
that while being a good approximation, is not exact; indeed, higher order terms in ( p!3D 
introduce odd powers of x — x'. For the GOP, the non-analytic behavior of ( p^ at the origin 
indicates that (|^ is not exact also in that case. 



Using the perturbative expression ([TsD, it is possible to express the rotation rate scaling 



factor V R directly as the RMS of the normalized conductance peak velocity 



where r^j = 1/4 and for the GOP the derivative is discrete. One can calculate R from the 
conductance peak data according to (^) and then use (|lTD or ([12D to determine D and thus 



the scaled parameter x that leads to universal correlations. A semiclassical calculation of R 



for a ballistic electron in a magnetic field {x = B) leads to an estimate of the correlation 
field = ^ (2mEA/h^y/\%/A), where A is the area of the dot and $o is the fiux 

unit. 

Finally, we computed the peak correlator for an asymmetric dot using (|I]), this time 
with r' 7^ F'^. The left and right leads each have their own channel correlation matrix with 
no correlation between them for a sufficiently large separation. Fig. ^ displays the GP 
universal predictions for cg for single- channel symmetric leads (F' = F^), and a comparison 
with Anderson model simulations. Eq. still provides a good fit but with ai = 0.37 ±0.04 
and 02 = 0.54 ± 0.04. We also find that Eq. holds but with h = 7/4 and 62 = 3 H. 
The rotation scaling factor is calculated from ([T5|) where now ri = 1/7 and r2 = 5/24 ||14|| . 
For asymmetric single-channel leads (F' 7^ F^'), cq is found to be intermediate between the 
symmetric single-channel leads correlator and the symmetric dot correlator. The latter is 
experimentally measurable in an asymmetric dot whose leads are made very asymmetric, 
since the conductance peak is then approximately proportional to the dominating width. 
For symmetric leads with many equivalent channels the peak correlator also approaches the 
width correlator | ]14| . 

In conclusion, we computed the universal conductance peak correlator and obtained a 
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good approximation of its functional form. We derived the rotation rate scaling which 
is useful when the level velocity scaling factor is not measurable and which can be directly 
calculated from the conductance peaks data using (|T5|) . Since the submission of this Letter, 
our prediction (|^) has been confirmed experimentally for broken time-reversal symmetry [jl6 



We also remark that the width correlator (B,!?!) is identical to the correlator of wavefunction 



intensity at a fixed spatial point which can be measured in microwave cavity experiments ||T7 
as a function of shape. This work was supported in part by the Department of Energy Grant 
DE-FG02-91ER40608. We acknowledge CM. Marcus and A.D. Stone for useful discussions. 
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FIGURES 

FIG. 1. Universal form of the width correlator (P) as a function of the scaled parameter (Q). 
Top: GP calculations (diamonds) and their best fit to (|^) (dashed). For comparison we also plot the 
GUP result on the left and the GOP result on the right (dotted). Insets: the universal ratio R/D 
computed with a finite Ax using simulations of the GP (diamonds with statistical errors) and the 
Anderson model (squares and circles). The dashed line in the left inset is the theoretical prediction 
([T^). Bottom: Anderson model simulations. Left: leads of 1 x 1 (pluses) and 4x4 (squares) points 
with a varying step potential. Right: a magnetic flux of is applied for 1x1 leads with an 
additional varying magnetic field (crosses) and 4x4 leads with a varying step potential (circles). 
Insets: typical width fluctuations for a single eigenfunction before the statistical averaging for the 
GP (solid) and Anderson model (dashed). 

FIG. 2. Universal form of the peak correlator in asymmetric dots with symmetric single-channel 
leads for both orthogonal and unitary symmetries. The GP predictions (dashed) are compared with 
Anderson model simulations (pluses and crosses as in Fig. 1). Shown by dotted lines is the width 
correlator which also describes the conductance correlator for highly asymmetric leads. 
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